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Abstract. Hadron physics represents the study of strongly interacting matter in all its manifes- 
tations and the understanding of its properties and interactions. The interest on this field has been 
revitalized by the discovery of new light hadrons, charmonium- and bottomonium-like states. In 
this talk I review the most recent experimental results from different experiments. 
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Introduction 

m . 

£> , The modern theory of the strong force, which binds quarks inside hadron matter, is quan- 

*sO ■ turn chromodynamics (QCD) fljj, |2fl. In the QCD framework only color-singlet states can 

exist due to confinement, and only some combinations of color states produce an attrac- 
tive potential, leading to a bound state. Therefore, the field of hadronic physics is the 
study of strong interaction in all visible matters and understanding of fundamental ques- 
tions in terms of QCD. Recently the experimental situation in this field appeared quite 
exciting. Especially, the situation changed dramatically in the last ten years. After the 
Belle Collaboration claimed the observation of a new and narrow resonance around 3.872 
GeV/c 2 , decaying into J/ipir + n + , named X(3872) [31], many similar charmonium and 
bottomonium states was found. Most of these states do not fit in the standard charmo- 
nium/bottomonium model. 



J-j ■ This review talk is devoted to the experimental signatures of these hadron spectroscopy, 

in particular, these exotic states. These results are from experiments at both e + e~ and 
hadron colliders. Most of the new charmonium states observed in the last few years 
are from B factory experiments BABAR and Belle. Many new observations of hadron 
spectroscopy also come from the r— charm factory experiments BESIII and CLEO-c, 
and from hadron colliders including experiments at Tevatron and LHC, as well as fixed 
target experiments, such as COMPASS. The analysis techniques adopted in the search 
for hadron spectroscopy depend on the production mechanisms at different experiments. 
In the B factories at e + e~ machine, the hadron can be produced in a B meson decay, 
e.g. B — > Kh, and in the s-channel direct production, e + e~ — > 7* — > h or with an 
initial state radiation (ISR) for hadron state carrying J PC = 1~~, and in the two photon 
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fusion production, as well as double quarkonium production, e.g. e + e~ — > (cc)(cc). At 
hadron collider, prompt productions coexist with the production in B meson decays. At 
a fixed-target experiment like COMPASS, the production mechanisms are mainly from 
three mechanisms: central production, diffractive dissociation and photoproduction. In 
this paper, for the reported experimental results, the first error and second error will be 
statistical and systematic, respectively, if they are not specified. 



2. Light meson decay and spectroscopy 

2.1 77 and 77' mesons and their excitation states 

The 77 and 77' mesons play an important role in understanding the low energy QCD. They 
are isoscalar members of the nonet of the lightest pseudoscalar mesons. Precision mea- 
surements on 77 and 77' would be very helpful and provide useful information in our un- 
derstanding of low energy QCD. 

The conversion decay of 77 — > e + e~e + e~ is important for the understanding of the 77 
coupling to the virtual photons and the calculation of the anomalous magnetic moment 
of the muon [4]. The KLOE Collaboration performed analysis based on data sample 
of 1.7 fb _1 collected at y/s = 1.02 GeV in the <f> meson mass region. By using the 
radiative decay <fi — > 777, the decay rate for the 77 — > e + e~e + e~(j) has been firstly 



obtained to be BR{ri ^ e + e~ e + e~ (7)) = (2.4 ±0.2 ±0.1) x 10" 5 1J. As a result, the 
measured branching ratio is fully radiation inclusive. The measurement is in agreement 
with theoretical predictions, which are in the range (2.41 — 2.67) x 10~ 5 [6-8]. FigureQ] 
shows the mass distribution of four electrons. 




Figure 1. M B 



mass distribution from data and fitting results. 



Precise measurements of the 77 and rf decay rates provide important information to test 
the Chiral Perturbation Theory (ChPT) [9]. The decays 77/77' — > 7r + 7r~7 are expected 
to get contributions from the box anomalies which proceed through a vector meson reso- 
nance, described by the Vector Meson Dominance ( VDM) [10]. According to the effective 
theory the 77/77' — > ir + Tr~-f processes are supposed to proceed both via a resonant con- 
tribution, mediated by the p meson, and a non-resonant direct term, connected to the box 
anomaly. Using a data sample corresponding to an integrated luminosity of 558 fb _1 at 



1 peak, the KLOE measured the ratio of 77 — > 7r + 7r 7 and 77 
> — > 777 [11]. The preliminary result is 



~7T° using the decay 



BTZ(rj — > 7r + 7r 7) 



= 0.1838 ±0.0005 ±0.0030. 



(1) 



BTZ{-q -> tt+tt-tt ) 
The result is the most precision measurement and in agreement with the recent results 
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from CLEO B12I1 . which differs by more than 3 a from the average of the previous re- 
sults fll . 

With new data sample of 225 million J/tp decay collected at the BESIII detector, 
the r\ and 77' decays can be studied via the charmonium decays into final states involv- 
ing 77/77' meson. The rare and forbidden decays of 77 and 77' could be reached at the 
BESIII experiment [14]. The BESIII Collaboration made a precision measurement of 
rf —> 3ns via the decays J/tp — > 77r + 7r~7r° and 77r 7r 7r° 11511 . The branching frac- 
tions are determined to be BTUyf ->■ ir+ir-ir°) = (3.83 ± 0.15 ± 0.39) x 10~ 3 and 



— ► 7T T 7T 7T 



II 



BK(r)' -> 7r 7r°7r°) = (3.56 ± 0.22 ± 0.34) x 10~ 3 . For the decay 77' 
the branching ratio is consistent with the CLEO-c measurement [16], and the precision is 
improved by a factor of four. For the decay rf — > 7r 7r°7r°, it is two times larger than the 
world average value IU3I1 . 
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Figure 2. The 7r + 7r~ (left) and 7r°7r° (right) invariant mass spectra with ty + tv~tt° 
(37r°s) in the 7/(1405) mass region. The solid curve is the result of the fit described 
in the text. The dotted curve is the /o(980) signal. The dashed curve denotes the 
background polynomial. 
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Figure 3. Results of the fit to (left) the /o(980)(7r + 7r~)7r and (right) 
/o(980)(7r°7r )7r° invariant mass spectra. The solid curve is the result of the fit de- 
scribed in the text. The dotted curve is the /i(1285)/??(1295) and 7/(1405) signal. The 
dashed curves denote the background polynomial. 



The spectrum of radial excitation states of isoscalar 77 and rf is still not well known. 
An important issue is about the nature of 77(1405) and 77(1475) states, which are not well 



established. By using the decays of J/ip 



77T"' 



and 77T 7r 7r° in BESIII, clear 



/o(980) signals are observed on both 7r + 7r~ and 7r°7r° spectra as shown in Fig. [2] the 
width of the observed /o(980) is much narrower (less than 10 MeV) than that in other 
processes 1 1311 . By taking events in the window of / O (980) on the 7r + 7r~ (7r°7r°) mass 
spectrum, evidence of h (1285)/?y(1295) is observed in the low mass region of /o(980)7r° 
as shown in Fig. [3] [15]. It is interesting that a clear peak around 1400 MeV is also 
observed on the mass of /o(980)7r° (see Fig. [3j. Angular analysis indicates that the peak 
on 1400 MeV is from 77(1405) -4 /o(980)tt decay. The BESIII Collaboration measured 
the product branching fraction of 7/(1405) production to be B7Z(J/ip — > 77/(1405)) x 
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/3"K(t7(1405) -» /o(980)vr ) x /3ft(/ (980) -> tt+tt-) = (1.50 ± 0.11 ± 0.11) x 10~ 5 
and BH(J/ip -> 77y(1405)) x /3ft(r?(1405) -> / (980)7r°) x Eft(/ (980) -> 7r°7r°) = 
(7.10 ± 0.82 ± 0.72) x 10~ 6 111, respectively. According to PDG values fH, one 
obtains the ratio £^(77(1405) ->■ f (980)ir°)/B1l(ri(UQ5) -► a (980)7r°) - 25%. It 
is the first time that we observe anomalously large isospin violation in the strong decay 
of 7y(1405) —> /o(980)7T° (for example, the ratio of the isospin violating rf —> 7r + 7r~7r° 
to the isospin covering 77' — > 7r + 7r~77 is about 0.8%). Following BESIII measurement, in 
reference [17!], the authors interpret this puzzle as an intermediate on-shell KK* + c.c. 
rescattering to the isospin violating /o(980)7r° by exchanging on-shell kaon. Further 
experimental study on the 7/(1405) parameters and identification of quantum number are 
needed at BESIII. 



2.2 Spin exotic light states 

COMPASS is a multi-purpose fixed-target experiment at the CERN Super Proton Syn- 
chrotron (SPS) aimed at studying the structure and spectrum of hadrons. One primary 
goal is the search for the spin-exotic mesons and glueballs. In a partial-wave analysis 



PC _ j-+ 



res- 



(PWA) of the pilot run data taken in 2004, a significant spin-exotic J 
onance was found at around 1660 MeV/c 2 in it~ii + tt~ final states produced in ir~ 
diffraction on a Pb target [18]. Its mass-dependent phase differences to the J PC = 2 h 
and 1 ++ waves are consistent with the highly debated 7Ti(1600) meson claimed in this 
channel by E852 and VES experiments [ 19, 20]. From a mass-dependent fit a resonance 
mass of (1660 ± W±° 64 ) MeV/c 2 and a width of (269 ± 21±<g) MeV/c 2 are deduced. 
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Figure 4. Top row. Intensities of major waves: 2 ++ l + [pir]D with (22(1320) (left), 
l ++ + [p7r]S with 01 (1260) (center) and 2~ + + [/ 2 7r]S with tt 2 (1670) (right). Bot- 
tom row: Intensity of the spin-exotic wave 1 h l + [p7r]P (left), and phase differences 
of this wave with respect to the l ++ + [pir]S (center), and the 2 h + [/27r]S' waves 
(right) (From j^). 



In 2008 COMPASS has acquired large data sets of diffractive dissociation of 190 GeV/c 
7r~ on a H2 target. In a partial-wave analysis (PWA) the isobar model 112 ill is used to de- 
compose the decay X~ —¥ ir~ir + ir~ into a chain of successive two-body decays. The 
spin-density matrix is determined by extended maximum likelihood fits performed in 20 
MeV/c 2 wide bins of the three-pion invariant mass. The intensity of the three domi- 
nant waves in the tt^tt+tt^ final state, l ++ 0+[p7r]5 , 2 ++ l+[p7r]Z), and 2 _+ + [/ 2 7r]S' 
, are shown in Fig. [4] top row. They contain resonant structures that correspond to the 
ai(1260), a2(1320), and 7^(1670), respectively 112211 . Of peculiar interest are the fit re- 
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suits for the spin-exotic wave in Fig. |4] bottom, left plot. The plot nicely illustrates the 
unprecedented statistical accuracy due to the large data set. The 1 v \ + \pix\P intensity 
features a broad bump, centered at 1.6 GeV/c 2 . In this mass region a rising phase with 
respect to the tail of the ai(1260) in the l ++ + [/97r]S' wave is seen (Fig. [4] bottom, cen- 
ter). As in Fig. 1 (bottom, right shows), the structure is phase locked with the 7^(1670) 
in the 2 h + [/27r]S' wave. This is consistent with the results obtained from a PWA of the 
pilot-run data taken with a Pb target 11811 . 

The CLEO-c Collaboration presents an analysis of V , (2<S') — > jXci 
cays in which they study the production of various r]^7T intermediate states [23]. They 
find evidence for an exotic rj'ir P-wave scattering amplitude at the level of 4 standard de- 
viations under a wide variety of model variations. The best fit to the data is obtained when 
the 7Ti7r amplitude is included. The ij'tt mass is described by a Breit-Wigner lineshape 
with a mass and width of 1670 ± 30 ± 20 MeV/c 2 and 240 ± 50 ± 60 MeV/c 2 , respectively, 
which is the first evidence of exotic state in charmonium decays. The result is consistent 
with that from other experiments [19, 20]. Now BESIII had collected 4 times of CLEO-c's 
sample at ijj{2S) peak B2411 . we expect to confirm CLEO-c measurement soon. 



^yrj^TT + TT de- 



2.3 New light mesons observed by BESIII 

In 2005, BESII observed a ?/7r + 7r~ resonance, X(1835), in the radiative decay J/i/j —> 
77/7r + 7r~ with a statistical significance of 7.7cr [25]. A fit to a Breit-Wigner function 
yielded a mass M = 1833. 7 ±6.1 ±2. 7 Me V/c 2 , a width T = 67.7 ±20.3 ±7.7 Me V/c 2 . 
The study was stimulated by searching for h glueball candidates which was predicted 
by the Lattice QCD [26]. With 225 million J/ip decay events collected by the BESIII 
detector, the X(1835) state has been confirmed with a statistical significance larger than 
20(T in the same analysis [27]. The mass and width are measured to be M = 1836.5 ± 
3.0±2:i MeV/c 2 and T = 190 ± 9+3* MeV/c 2 as shown in shown in Fig. [5] The mass 
of the X(1835) is consistent with the BESII result, but the width is significantly larger. 
A simple angular analysis indicates that quantum number of X(1835) is consistent with 
a pseudoscalar assignment, but the others are not excluded. Meanwhile, two resonances, 




1.4 1.6 l.S 2.0 2.2 2.4 2.6 2.8 
M(7t*jnf)(GeVA?) 



Figure 5. (a) The rj'-K + iT^ invariant-mass distribution, (b) mass spectrum fitting 
with four resonances, here, the dash-dot line is contributions of non-77' events and the 
7/ tv + -k~ 7T° background and the dash line is contributions of the total background and 
non-resonant rj'-K + -K~ process. 
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the X(2120) and the X(2370) are observed with statistical significances larger than 1 2a 
and 6.4cr, respectively, as shown in Fig. [5] The masses and widths are measured to be: 



X(2120) 



M = 2122.4 ± 6.7i^7 MeV/c 2 , r = 83 ± l&t\\ MeV/c 2 



(2) 



X(2370) 



M = 2376.3 ± 8.7+^3 MeV/c 2 , r = 83 ± 17+g MeV/ 



(3) 



In the mass spectrum fitting in Fig. |5jb), possible interferences among different reso- 
nances and the non- resonant process are not taken into account. To determine the spin 
and parity of the X(1835), X(2120) and X(2370), and to measure their masses and 
widths more precisely, a PWA must be performed, which will be possible with the much 
higher statistics J/ip data samples planned for future runs of the BESIII experiment. 
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Figure 6. Invariant-mass distributions for the selected events: left is the invariant-mass 
spectrum of t\-k\ right: Results of the fit to the r;n + Tv~ mass distribution for events with 
either r]-K + or tj-k~ is in the ao(980) mass window. 



The BESIII experiment also reports an analysis of J ftp — > u)r\-K + ix~ . A structure 
around 1.8-1.9 GeV/c 2 in the i]7r + ir~ mass spectrum is observed [28]. In the analysis, 
a clear ao(980) signal is seen in the r/ir mass spectrum as shown in Fig. [6] (left). The 
?/7r + 7r~ mass spectrum for events where either M(rjir + ) or M(rjir^ ) is in a 100 MeV/c 2 
mass window centered on the ao(980) mass is shown in Fig. [6] (right). Three peaks are 
observed on the rjTT + ir^ mass spectrum, two of them in the low mass side are /i(1285) 
and ?/(1405). The structure near 1.8 GeV/c 2 is the first observation. A fit to the three 
signal peaks leads to M = 1877.3 ± 6.3±^ MeV/c 2 and T = 57 ± 12±f MeV/c 2 
for the AT(1870) structure with the statistical significance of 7.2o\ Here the three signal 
peaks are parametrized by Breit-Wigner functions convolved with a Gaussian resolution 
function and multiplied by an efficiency curve, which are both determined from signal 
MC samples and fixed in the fit. Whether the resonant structure of X(1870) is due to the 
AT(1835), the 772(1870), an interference of both, or a new resonance still needs further 
study such as a PWA that will be possible with the larger J/ip data sample. 

2.4 Proton-anti-proton mass threshold enhancement 

An strong pp mass threshold enhancement was first observed by the BESII experiment in 
the decay process J/tp — > ^pp [29] and was confirmed by the CLEO-c experiment [30]. 
With the BESIII data, a PWA analysis is performed to determine the parameters of the 
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pp mass threshold structure, which we denote as X(pp) Bill . Only events of M p p < 2.2 
GeV/c 2 are considered in the PWA. Four components, the X(pp), / 2 (1910), / (2100) 
and ++ phase space (PS) are included in the PWA fit. The intermediate resonances are 
described by Breit-Wigner propagators, and the parameters of the /2(1910) and Jo (2100) 
are fixed at PDG values. In the fit, the pp final state interaction (FSI) effect is also consid- 
ered by using the Julich formulation [32]. Figure Q shows comparisons of the mass and 
angular distributions between the data and the PWA fit projections. In the optimal PWA fit, 
the X(pp) is assigned to be a h state, which is 6.8rr better than other J PC assignment. 
In the fit, a Breit-Wigner and S-wave final state interaction (1=0) factor can well describe 
the pp mass threshold structure. The PWA fits are also performed without the correction 
for FSI effect. The corresponding log-likelihood value worsen by 25.6 than those with FSI 
effect included. With the inclusion of Julich-FSI effects, the mass, width and product BR 
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Figure 7. Comparisons between data and PWA fit projection: (a) the pp invariant mass; 
(b)-(d) the polar angle 8-, of the radiative photon in the J/tp center of mass system, the 
polar angle 6 V and the azimuthal angle 4> p of the proton in the pp center of mass system 
with Mpp — 2 Alp < 50 MeV/c 2 , respectively. Here, the black dots with error bars 
are data, the solid histograms show the PWA total projection, and the dashed , dotted , 
dash-dotted and dash-dot-dotted lines show the contributions of the X(pp), ++ phase 
space, /o(2100) and /2(1910), respectively. 



for the X(pp) are measured to be: M = 1832±,t 9 (siat.)±^(syst.) ± 19 (model) MeV/c 2 , 
T = 13±39(stat.)±}°(syst.)±4(model)MeV/c 2 (a total width of T < 76 MeV/c 2 at 90% 
C.L.) and BH(J/ip -» 1 X)B1l{X -> pp) = (9.0±$;^ (staL)±^J (syst) ± 2.3(model)) x 
10~ 5 , respectively, where the third error are uncertainty due to choosing different model 
of FSI effect JUH. 

The tp(2S) — > jpp decay is also studied by using 106 million tp(2S) decay events 
collected by the BESIII detector [31]. The pp mass spectrum of the surviving events 
is shown in Fig. Oa). Besides the well known r/ c and Xcj peaks, there is also a pp 
mass threshold excess relative to phase space. However, here the line shape of the 
mass spectrum in the threshold region appears to be less pronounced than that in J/tp 
decays. A PWA on the selected tp(2S) — » jpp which is similar to that applied for 
J/tp — > rypp decay was performed to check the contribution of X(pp) in ip(2S) de- 
cays and to measure the production ratio between J/i/j and tp(2S) radiative decays, 
R = BTZ{ip{2S) -> jX(pp))/BK(J/ip -> ~/X{pp)). Due to limited statistics of^{2S) 
events, in the PWA, the mass and width of X(pp) as well as its J PC were fixed to the 
results obtained from J/tp decays. Figure [8] (b) shows comparisons between data and 
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Figure 8. (a) The pp invariant mass spectrum for the selected ip(2S) — > "/pp candi- 
date events; the open histogram is data and the dashed line is from a ip(2S) — > jpp 
phase-space MC events(with arbitrary normalization), (b) Comparisons between data 
and PWA fit projection for pp mass spectrum, the representations of the error bars and 
histograms are same as those in Fig. [7] 



MC projections for the pp mass spectrum. The determined product B7Z and R value are 
BTl{ip{2S) -> jX)BK{X -> pp) = (4.57±0.36(stat.)tJ; U 7(syst.) ± 1.28(model)) x 
10~ 6 and R = (5.08lJ];^(stat.)i°g(syst.) ± 0. 12 (model) )%, respectively. It is sup- 
pressed compared with the so called 12% rule. 



3. Charmonium spectroscopy and decays 

3.1 T] C (1S) and r/ c (2S) 

3.1.1 rj c (lS) resonance via ip(2S) —> "frj c decay 

Precise measurement of Ml transition of i/j(2S) is important for us to understand the 
QCD in the relativistic and nonperturbative regimes. The i/j(2S) — * ^r\ c transition is also 
a source of information on the rj c mass and width. There is currently a 3.3<r inconsistency 
in previous r\ c mass measurements from J/ip and tl>(2S) — > ^r\ c (averaging 2977.3 ±1.3 
MeV/c 2 ) compared to 77 or pp production (averaging 2982.6 ±1.0 MeV/c 2 )[13]. The 
width measurements also spread from 15 to 30 MeV, it is around 10 MeV in the earlier 
days of experiments using J/ip radiative transition II35U36I1 . while the recent experiments, 
includi ng p hoton-photon fusion and B meson decays, gave higher mass and much wider 
width l37l - J40ll . The recent study by the CLEO-c experiment [41], using both ip(2S) 
and J/ip — > jr/ c decays, and pointed out there was a distortion of the r\ c line shape. 
The CLEO-c attributed the r\ c line-shape distortion to the energy-dependence of the Ml 
transition matrix element. In the J/ip —> ^r\ c from CLEO-c, the distorted r\ c lineshape 
can be described by the relativistic Breit-Wigner (BW) distribution modified by a factor 
of E^ together with a dumping factor to suppress the tail on the higher photon energies. 
The KEDR Collaboration did the same thing but tried different dumping factor 04211 . 

Based on the data sample of 106 M ip(2S) events collected with the BESIII detector, 
the rj c mass and width are measured from the radiative transition i/j(2S) —> -fr] c 143I1 . 
The j] c candidates are reconstructed from six exclusive decay modes: K + K~ir + ir~ir°, 
K + K~ir°,K s K + ir~ir + Tr~ , K s Kir, i]ir + ir~ , and 3(7r + 7r~), where K s is reconstructed 
in 7r + 7r~ mode, r\ and tt° from 77 final states. For a hindered Ml transition the matrix 
element acquires terms proportional to E^, which, when combined with the usual E^ 
term accounting for the wave function mismatch between the radial excited tp(2S) and 
the ground-state ^(lS") transitions, lead to contributions in the radiative width propor- 
tional to E 1 . Thus, the r\ c lineshape is described by a BW modified by EJ. convoluted 
with a resolution function. It is important to point out that the interference between r\ c 
and non-resonance in the signal region is also considered. The statistical significance of 
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the interference is 15o\ This affects the ?y c resonant parameters significantly. By assum- 
ing all non-resonant events interfere with the rj c , the BESIII analysis obtains r\ c mass and 
width, M = 2984.3 ± 0.6 ± 0.6 MeV/c 2 and T = 32.0 ± 1.2 ± 1.0 MeV, respectively. 
Two solutions of relative phase are found for each decay mode, one represents construc- 
tive interference, the other for destructive. Regardless which solution one takes, the mass, 
width of the r/ c and the overall fit are always unchanged. The relative phases for con- 
structive interference or destructive interference from each mode are consistent with each 
other within 3a, which may suggest a common phase in all the modes under study. The 
fit with common phase for each mode shows the relative phase <p = 2.40 ± 0.07 ± 0.08 
rad (constructive) or <\> = 4.19 ± 0.03 ± 0.09 rad (destructive). The physics that the inter- 
ference phase is the same for the six channels is yet to be understood. Figure |9]shows the 
fit results in the six r\ c decay modes. 
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Figure 9. The invariant mass distributions for the decays 

K + K-TT + n-n ,K + K-TT ,K 3 K + TT-n + TT~,K s K-K,rin + TT- and 3(7r + 0, 
respectively. Solid curves show the fitting results; the fitting components (rj c sig- 
nal/non-resonance/interference) are shown as (dashed/long-dashed/dotted) curves. 
Points with error bar are data, shaded histograms are (in green/yellow/magenta) for 
(continuum/other 77c decays/other tp(2S) decays) backgrounds. 



With precise measurement of the r\ c mass, one can obtain the hyperfine splitting, AMhj(lS) c 
M(J/ip)—M(f] c ) = 112.6±0.8MeV, which agrees with the quark model prediction h4l , 
and will be helpful for understanding the spin-dependent interactions in hidden charmo- 
nium states. 



3.1.2 Observation ofip(2S) -J> 7^(25) 

The first radially excited S-wave spin singlet state in the charmonium system, r/ c (2S), was 
observed by the Belle experiment in the decay process B ± — > K ± r] c (2S), ri c (2S) — > 
K s K ± tt t [45]. It was confirmed by the CLEO 10] and BABAR Q experiments 
in the two-photon fusion process e + e~ — > e + e~ (77) , 77 — > rj c (2S) — > K S K^^ 
and by the BABAR experiment in the double-charmonium production process e + e~ — > 
J /4>(cc) 114811 . The only evidence for rj c (2S) in the ip(2S) — > jr/ c (2S) decay was from 
Crystal Ball Collaboration Il49ll by looking at the radiative photon spectrum. Recently, 
the CLEO-c experiment searched for the tp(2S) — > 7?/ c (2S f ) signal with i] c (2S) exclusive 
decay into 1 1 modes by using 25.9 M tp' events, and no evidence found. Product branch- 
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ing fraction upper limits are determined as a function of T(rj c (2S)) for the 1 1 individual 
modes |50]. 




»fiL<P«v*?> 



Figure 10. Preliminary fitting of the mass spectrum for rj c (2S) — > K s K 7r T . 



The BESIII experiment searched for the Ml transition i)(2S) — > -fi] c (2S) through the 
hadronic final states K s K ± tt t . A bump is observed around 3635 MeV/c 2 on the mass 
spectrum as shown in Fig.[l0] In order to determine the background and mass resolution 
using data, the mass spectrum range is enlarged (3.47 ~ 3.72 GeV/c 2 ) to include \ci 
and Xc2 events. The resonances Xci and \ci are described by the corresponding Monte 
Carlo (MC) shape convolved a Gaussian which takes account the small difference on the 
mass shift and resolution between data and MC. So the mass resolution for the rj c {2S) 
in the fitting is fixed to the linear extrapolation of the mass resolutions from the Xci and 
Xc2 signals in data. The lineshape for n c (2S) produced by such the Ml transition is 
described by (E?, x BW(m) x damping (-Ely)) <S> Gauss(0, a) where m is the invariant 

mass of KgK^nT , E 1 — *' a — is the energy of the transition photon in the rest frame 

of ip(2S), damping (E 7 ) is the function to damp the diverging tail raised by E^ and 
Gauss(0,a) is the Gaussian function describing the detector resolution. The possible 
form of the damping function is somewhat arbitrary, and one suitable function used by 
KEDR for a similar process is (42 ] 



El 



E 7 E + (EjE — Eq) 



where Eq 



2 m 



i>(2S) 



is the peaking energy of the transition photon. In the fit, 



the width of n c (2S) is fixed to PDG value. From the fit to the data, a signal with 
a statistical significance of 6.5 standard deviation is observed which is the first obser- 
vation of the Ml transition ip(2S) — > jrj c (2S). The measured mass for i] c (2S) is 
3638.5 ± 2.3 ± 1.0 MeV/c 2 . The measured branching ratio is BK(t/j{2S) ->■ -/r] c {2S)) x 

BK(r] c (2S) ->• KsK^w^l = ( 2 ' 98 ± °- 57 ± °' 48 ) x 10 ~ 6 - To S ether with the BABAR 
result BK{r] c (2S) -> KKir) = (1.9 ± 0.4 ± 1.1)% JH, the Ml transition rate for 
ip(2S) -> JVc(2S) is derived as BTl{ip{2S) ->• 7r? c (25)) = (4.7 ± 0.9 ± 3.0) x 10~ 4 , 
which is in agreement with the prediction of potential model calculations |52]. 

Recently, the Belle experiment measured rj c (lS) and i] c (2S) resonant parameters in the 
decay of B ± — > K ± (K s Ktt)° by considering the interference between i] c (lS)/i] c (2S) 
decay and non-resonant in the S ± decay 15311 . Meanwhile, the BABAR experiment also 
updated the analysis of e + e~ — > e + e _ (77),77 — > r) c (lS)/i] c (2S) — > (K s Ktt)° and 
K + K~ir + ir~Tr° modes [54]. Table [T] shows the summary of the typical production pro- 
cesses (ijj(2S) radiative decays, 77 fusion and B decays) for rj c (lS)/r] c (2S) and corre- 
sponding resonant parameter measurements. These results indicate that the rj c (lS) pa- 
rameters agree well from different production processes. 
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Table 1. Comparison of the mass and width for r/ c (lS)/rj c (2S) in different production 
processes, ip' — >■ "/T]c(lS)/r] c (2S), B ± — > K ± r/ c (lS)/'q c (2S) and 77 fusion, from 
different experiments. The PDG values are only world average from earlier results. 
. For the time being, the most precise measurements for r/ c (lS) resonance are from 
BESIII, while these for r/ c (2S) resonant parameters are from BABAR in 77 fusion. 





BESIII 
4>' -* fVe/Vc(?S) 


Belle 153] 
B decays 


BABAR [54J 
77 fusion 


PDG 2010 [13] 


M(r) c (lS)) MeV/c 2 


2984.3 ±0.6 ±0.6 


2985.4 ± 1.5^2 


2982.2 ± 0.4 ± 1.4 


2980.3 ± 1.2 


r(jj c (lS))MeV 


32.0 ± 1.2 ± 1.0 


35.1 ±3.1 + i° 


32.1 ± 1.1 ± 1.3 


28.6 ± 2.2 


M(i7 c (2S))MeV/c 2 


3638.5 ± 2.3 ± 1.0 


3636.1+.? lt° 2 i 


3638.5 ±1.5 ±0.8 


3637 ± 4 


r(jj c (2S))MeV 


12 (fixed) 


fi R +8.4 + 2.6 
D ' D -5. 1-0.9 


13.4 ±4.6 ±3.2 


14 ±7 



3.2 New results on h c from BESIII 

The BESIII Collaboration reported the results on the production and decay of the h c using 
106M of ■0(25') decay events in 2010 [55], where they studied the distributions of mass 
recoiling against a detected ir° to measure ip(2S) — > ir°h c both inclusively (El-untagged) 
and in events tagged as h c — > jric (El -tagged) by detection of the El transition photon. 
In 201 1, 16 specific decay modes of r\ c are used to reconstruct r\ c candidates in the decay 
mode of h c — > ^r\ c . Figure [TTI (left) shows the ir° recoiling mass for the sum of the 16 
r\ c decay modes. Fits to the 16 tt° recoil-mass spectra are performed simultaneously that 
yields M(h c ) = 3525.31±0.11±0.15MeV/c 2 and T(h c ) = 0.70 ±0.28 ±0.25 Me V/c 2 , 
respectively. These preliminary results are consistent with the previous BESIII inclusive 
results and CLEO-c exclusive results. 




18 n 3.50 3.52 354 -3.5 

7i recoil mass (GeVAr) 



Figure 11. Left: The tv° recoiling mass for the sum of 16 i] c decay modes and fitting 
results. Right: The invariant mass distributions for the sum of 16 r\ c decay modes, and 
fitting results. All are preliminary results. 



The centroid of the 3 Pj states (Xc05Xci,Xc2) is known to be (M( 3 Pj)) = [5M ( 3 P 2 ) + 
3M ( 3 Pi) + M( 3 P )} = 3525.30 ± 0.04 MeV HI. If the 3 Pj states centroid mass 
(M( 3 Pj)) is identified as the mass of M( 3 P), then BESIII observes the hyperfine split- 
ting as AM hf (lP) cc = -0.01 ± 0.11 ± 0.14 MeV which agrees with zero. 

The BESIII Collaboration also looked at the 7] C (1S) mass distributions in the exclusive 
h c —> jr) c (lS) decay modes. Figure [TT] (right) shows the distribution of invariant mass 
distribution from exclusive hadronic decays of i] c (lS), the shape of the r\ c is symmetric in 
h c (lP) — > jric(lS) radiative transition and no distortion observed as in J/i/j/i/j' radiative 
decays 114 ll 14311 . A maximum likelihood fit is perform to extract the rj c (lS) parameters. 
In the fit the signal shape is described by a BW (corrected by E 3 ) convoluted with a 
Gaussian resolution function. The preliminary results for the r] c (lS) parameters are M = 
2983.6 ±1.1 MeV/c 2 and T = 36.6 ± 3.15 Me V/c 2 , where the errors are statistical error 
only. The results are in agreement with those from fj)(2S) — > jri c (lS) analysis. 
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3.3 Determination of ip(3770) parameters at KEDR 



A measurement of the 0(3770) meson parameters has been reported by the KEDR experi- 
ment based on data samples during the scans of the center-of-mass energy range from 3.67 
to 3.92 GeV at the VEPP-4M e+e - collider in 2004 and 2006 S. The observed mul- 
tihadron cross sections were fitted as a function of the center-of-mass energy using some 
assumptions about the behaviour of then non-resonant form factor. The data points cor- 
rected for the detector efficiency together with several fitting curves are shown in Fig. [? ] 
(left). Unlike Ref. [57], KEDR experiment did not observe any shape anomaly. The inter- 
ference of resonant and nonresonant production essential in the near-threshold region has 
been taken into account as suggested in Ref. [58]. The analysis assumes the domination 
of ij)(2S) in the non-0(377O) part of the D-meson form factor. It is also worth noting 
that correct accounting of the resonance-continuum interference can help in solving the 
0(3770) non-DD puzzle Q. 





3720 3740 3760 3780 3800 3820 3840 3860 



Figure 12. Left: Cross section of e + e~ — >• hadrons vs. cm. energy in the vicinity of 
0(3770) with the light quark, r and QED backgrounds subtracted. The curves are the 
result of a simultaneous fit for different models. Right: The same plot as left one, solid 
and short-dashed curves correspond to two VDM solutions. Resonant and non-resonant 
parts are presented separately. 



It is known that for two interfering resonances the ambiguity can appear in the reso- 
nance amplitudes and the interference phase [60]. Figure [12] (right) shows the results of 
the two solutions corresponding to vector dominance models (VDM). The KEDR exper- 
iment presented the results on the mass ad width of '0(3770) are M = 3779. ll^ ± 



0.618:1 MeV/c 2 and T 



25.2tti±0.5; 



I82 MeV, where average values of the mass and 



total width of the two VDM solutions are used and the third error arises from the model 



dependence. The results agree with those from BABAR [61] also taking into account 
interference and disagree with all results obtained ignoring this effect including that by 
BESII 1620 . The KEDR also obtained two possible solutions for the 0(3770) electron 
partial width 



1. r„ 



U7t% ± 13±i„ eV, 



2. T ee = 415+12 ± 38+iS eV 



58 



-10 



Most of the potential models strongly support the first solution and can barely tolerate the 
second one [63-65]. 
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Table 2. Comparison of the X(3872) mass measurements from different experiments, 
and world average. 



Experiment 


X(3872) massMeV/c 2 


CDF 


3871.61 ±0.16 ±0.19 


BABAR(B±) 


3871.4 ±0.6 ±0.1 


BABAR(S°) 


3868.7 ±1.5 ±0.4 


DO 


3871.8 ±3.1 ±3.0 


Belle (full dataset) 


3871.84 ±0.27 ±0.19 


LHCb 


3871.96 ±0.46 ±0.10 


World Average 


3871.67 ±0.17 



3.4 News on charmonium-like states 

3.4.1 X(3872) 

The X(3872) resonance is the forerunner of the new charmonium family. It was dis- 
covered by Belle [3] in the B — > KX decay, with X — > J/jjrrr + Tr~, and confirmed by 
BABAR S@l, CDF H and DO [69}, and now by CMS Sill and LHCb [72]. The 
quantum number J PC of the X(3872) state is not yet fully established. A study of the 
angular decay distributions performed by the CDF experiment gave nonnegligible prob- 
abilities only for 1 ++ and 2 h [68], where C = + is also confirmed by the observation 
of the X — > J J -07 decay M73I1 . Table [2] summarizes the measurements of the X(3872) 
mass. The new world average is M — 3871. 67 ± 0.17. The close proximity of the world- 
averaged mass to the D*°D° (m D ,a + m D o = 3871.79 ± 0.30 MeV/c 2 ) mass threshold 
has engendered speculation that the X(3872) might be a loosely bound D*°D° molecule 
or a tetraquark state [74]. Despite a large experimental effort, the nature of this new state 
is still uncertain and several models have been proposed to describe it. 

With full dataset (711 fb _1 ) collected by the Belle detector, the difference in masses 
of the X(3872) states produced in B+ ->• K+n+TT-J/tp and B° ->■ K°-k + tt- J/ip 
decays is determined to be J75ll: AM X(3872) = (-0.69 ±0.97 ±0.19) Me V/c 2 , which is 
consistent with zero and disagrees with theoretical predictions based on a diquark model 
for the X(3872) B76I1 . We conclude from this that the same particle is produced in the 
two processes and use a fit to the combined neutral and charged B meson data samples 
to determine: M x(3872 ) = (3871.84 ± 0.27 ± 0.19) MeV/c 2 . This result agrees with the 
current world-average value of 3871.67 ± 0.17 MeV/c 2 as listed in Table|2] The width of 
the X(3872) signal peak is consistent with the experimental mass resolution and we set 
a 90% CL limit on its natural width of Tx(3872) < 1-2 MeV, improving on the previous 
limit of 2.3 MeV. The details of the analysis can be found in [75]. The production ratio of 
X(3872) in the B° and B + meson decays are measured to be : 

= g -> gy (3872)) = 
v ; B(B+ -> K+X(3872)) ' v ; 

This value is above the range preferred by some molecular models for the X(3872): 
0.06 < R{X) < 0.29 & The BaBar result for this ratio is R(X) = 0.41 ± 0.24 ± 
0.05 B67I1 . Both BABAR 1 78] and Belle [75] experiments had searched the charged partner 
of the X(3872) decaying into 7r + 7r° J/tp by using B — > Ktt + tt°J/iP decay, and shows no 
evidence for a charged partner to the X(3872) decaying as X + — > p + J/ip. The current 
data favor isospin 7 = 0. However, the close proximity of the D* a D° threshold may 
induce large isospin violations, as pointed out by N. A. Tornqvist [79]. Unfortunately, 
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in Belle experiment, an attempt to discriminate between J PC = 1 ++ and 2 h with an 
angular analysis of X(3872) — > ir + ir~ J ftp didn't brought to any conclusion, due to the 
lack of statistics [75]. 

Radiative decays of the X(3872) are important in understanding its nature. The de- 
cay of X(3872) -»■ J/V>7 established its charge parity to be +1. Both BaBar Hi and 
Belle [73] searched for the the X(3872) -> J/V> 7 and X(3872) -> ip(2S)j decays, 
that are predicted to dominate for a molecule. The BABAR Collaboration [80] show that 
B7e(X(3872) -► ip(2S)j) is almost three times that of BK(X(3872) -> J/ip-y), while, 
the Belle Collaboration find no evidence for X(3872) — > ip(2S)^, and determine the ratio 



to be : R 



6(X(3872)^t// 7 ) 



< 2.1 (at 90% C.L.) [73]. The possibility of a tp(2S)-y 
The X(3872) state may not have a large cc admix- 



B(X(3872)-*.J/V7) 

dominance seems to be excluded 

ture with a D*°D° molecular component as was expected on the basis of the BABAR 

result H. 

First observation of the X(3872) state at the LHC were reported by both CMS J70[|7lll 
and LHCb [72] experiments in the decay of J/ipTT + n~, using data collected in 2010. With 
an integrated luminosity of 34.7 pb _1 collected by the LHCb experiment, the production 
of the X(3872) particle is observed in pp collisions at y/s = 7 TeV. Figure [TH (left) 
shows the results. The masses of both the X(3872) and ip(2S) are measured to be : 
™*(3872) = 3871.95 ± 0.48 ± 0.12 MeV/c 2 and m^ (2S) = 3686.12 ± 0.06 ± 0.10 
MeV/c 2 . The CMS Collaboration performed a measurement of the X(3872) with data 
set collected in 2010 and 201 1, and found about 5 3000 candidates in 896 pb _1 , as shown 
inFig.[[3](right). 



-HCb i 

1=7 TeV 




tf ^^^ W ^ 



3600 3700 



M(J/>ji*ir}[MeWc*] 



i 




mlJ/'l' * x'X] GeWc* 



Figure 13. Left from LHCb: Invariant mass distribution of J /iJ)-k + -k~ (black points 
with statistical error bars) and same-sign J/ipir ir (blue filled histogram) candidates. 
The solid red curve is the result of the fit. Right from CMS: Invariant mass distribution 
of J/ip-K + -K~ from. The inset shows a zoom of the X(3872) region. 



3.4.2 1 family 

In the family of the charmonium-like states, the 1~ _ states are among the most easi- 
est to be found at e + e~ collider, since they can be produced through the ISR mecha- 
nism. The F(4260) is the first one found by BABAR Collaboration §£M. Now four 
states were discovered at the B-factories JUrH: the F(4008) and the F(4260) de- 
caying to J/^tt+tt-, the F(4350) and the y(4660) decaying toj/>(2S)7r + 7r-. The 
y(4260) has been also confirmed by the CLEO-c experiment F85L 18611 that could exten- 
sively produce it by running at open charm energy region. This allowed to observe also 
the Y — > J/ilm°ir° and the Y — > J/ijjK + K~ modes. The Belle experiment performs 
an analysis to compare the F(4260) —> J/i/jtt + tt^ and the F(4260) —> J/i/jt^ ^ rates. 
According to the isospin symmetry, the second one should be half of the rate of the first 
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J/ , 07r O 7r as function of mass have been measured, 
and by fitting to the cross-sections data, they obtain T ee x £>7£(F(4260) — > J/-07r°7r°) = 
(3.19lti'53io 35) e ^' which is consistent with the isospin expectation by comparing to the 

measured T ee x BTZ(Y{A260) -> J/V-tt+tt") d. 

3.4.3 Y (4140) atLHCb 

The CDF experiment reported a 3. 8cr evidence forthe Y (4140) state in the B — > KY(Y — > 
J/ip(j>) decay using pp data collected at the Tevatron (y/s = 1.96 TeV) [87]. A prelim- 
inary update of the CDF analysis with 6.0 fb _1 leads to observation of F(4140) with 
with a statistical significance of more than 5a 113811 . The mass and width were measured 
to be 4143^30 ± °- 6 MeV and 15 -6°i 4 ± 2 - 5 MeV/c 2 . The relative decay rate was 
determined to be B1Z{B+ -> YK+) x BK(Y -> J/^(j})/BTZ{B+ -> J/ip<t>K+) = 
0.149 ± 0.039 ± 0.024. However, the LHCb Collaboration perfrom the most sensitive 
search for the narrow Y"(4140) — > J/ip(f> state in _B + — > J/ip<f> decays by using 0.37 
fb _1 data. They did not confirm the existence of such a state. An upper limit on the 
Bn Ut -+ YK+ ) x BTZ{Y -» J/^(j>)/BTZ{B+ -> J/i;(f>K + ) < 0.07 at 90% C.L. is 
set (89J. The result disagrees at the 2.4er level with the CDF measurement. 

4. Bottomonium-like states 

The Belle experiment has collected a large sample of e + e~ collision at the center-of- 
mass energy near the Y(55) resonance, which lies above the B S B S production threshold. 
Many unexpected non-B s B s decays of Y(55) have been observed. In particular, anoma- 
lously large rates for dipion transitions T(55) — > T(nS)Tr + ir~ (n = 1, 2,3) have been 
observed [90]. Assuming these signals are attributed entirely to the T(55) decays, the 
measured partial decay widths T(T(55) — > T(nS)n + Tr~) ~ 0.5 MeV are about two 
order of magnitude larger than typical widths for dipion transitions of T(25), T(35) and 
T(45). 

Recently the CLEO-c experiment observed the process e + e~ — > h c (lP)ir + Tr~ at a 
rate comparable to the process e + e~ — > J/tpir + Tr~ by using data sample taken near 
y/s = 4170 MeV and found evidence of an even high transition rate at the F(4260) 
energy [91]. This implies that the ht,(nP) production might be enhanced in the region of 
the T(55), which may exist an exotic resonance Yj, analogue of the F(4260). 

Using the full Y(55) data sample with the integrated luminosity of 121.1 fb _1 col- 
lected near the peak of Y(55) with the Belle detector, they observe the hb(lP) and 
hb(2P) in the missing mass spectrum of tt + tt~ pairs. The ir + Tr~ missing mass is de- 



fined as MM(tt + tt ) = J (E c . m . — E* +7r _) 2 — p* 2 ^-, where E c _ m , is the center-of- 

mass (cm.) energy, E* + _ and p* , _ are the ir + ir~ energy and momentum measured 
in cm. frame. The details of the analysis can be found in [92, 93]. The MA1(tt + it~) 
spectrum with the combinatorial background and K s contributions subtracted, and the 
signal function resulting from the fit overlaid, are shown in Fig. [14] The significance 
of the hf,(lP) and hb(2P) signals which includes the systematic uncertainty is 5.5c 
and 11.2cr, respectively. This is the first observation of the h\,(\P) and hb(2P) spin- 
singlet bottomonium states in the process e + e~ — > hb(nP)n + ir~ at the Y(55) en- 
ergy. The measured masses and cross-sections relative to the e + e~ — > T(2S)tt + it~ 
cross-section are M = 9898.25 ± 1.06 j^g? MeV/c 2 , R = 0.407 ± 0.079±g;g^ for 
the h b (lP) and M = 10259.76 ± 0.64±^;g| MeV/c 2 , R = 0.78 ± 0.09±g;fg for the 
hb(2P). The mass of the spin-singlet state is consistent with the center-of-gravity of the 
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10.2 10.4 

MM(jiV), GeV/c 2 



Figure 14. The MM(-k + -k~) spectrum with the combinatorial background and K 8 
contributions subtracted (dots with error bars) and signal component of the fit function 
(solid histogram). The vertical dashed lines indicate the boundaries of the fit regions. 



corresponding \bj states. The hyperfine splitting is AM (IP) h f — 1.62 ± 1.52 MeV/c 2 
(AM{2P) HF = QA8,±Hl MeV/c 2 ). 

The large R values indicate that an exotic state, Y&, may be around T(5S) energy, 
and the exotic decay violates the suppression of heavy quark spin-flip. Belle experi- 
ment studied the resonant substructure of T(5S) — > hb(nP)ir + ir~ decays [92]. Since 
the low statistics and high background, a Dalitz plot analysis is impossible. Belle stud- 
ied performed the one-dimensional distribution of M(hb(nP)ir), where M(hi,(nP)ii + ) 
(M(hb(nP)%~)) is defined as a missing mass of the opposite-sign pion, M M (ir~) (MM(w + )). 
The hb(nP) signal yields are measured as a function of the MM(7r ± ) by fitting the 
MM(it + tt~) spectra in the bins of A/Af(7r ± ). Figure [15] (a) and (b) show results of 
the fits for the h^llP) and hi,(2P) yields as a function of MM (tt). A clear two-peak 
structure is seen without any significant non-resonant contribution. By assuming that the 
spin-parity for both structures is J p = 1 + , a \ 2 fit to the MM(tt) distributions is per- 
formed. In the fit, two P-wave BW amplitudes and a non-resonant contribution is used. 
The results of the fit are shown in Fig. Q3] and are summarized in Table [3] The non- 
resonant amplitude is found to be consistent with zero. In the fit, the hypothesis of two 
resonances is favored over the hypothesis of a single resonance (no resonances) at the 
7.4 a (17.9 a) level. Therefore, the two charged structures are named as Z;,(10610) and 
•Z&(10650), respectively. The parameters of the Zb(10610) and Z;,(10650) obtained in 
the fit of hi,(lP) and hb(2P) are consistent with each other. 




> 1 7500 

i 

o 15000 




1 




2 12500 






J 




iS 10000 


- 


'It 






7500 


- 


l l 






5000 


T 


1 


irl 


2500 


rl 




It 





" , 1 




h 



10.6 10.7 

MM(n), GeV/c 2 



10.6 10.7 

MM(n), GeV/c 2 



Figure 15. Left: the yield of the ht(lp) as a function of MM(ji) (points with error 
bars) and results of the fit (histogram). Right: the yield of the hb(2P) as a function of 
MM (it) (points with error bars) and results of the fit (histogram). 
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Table 3. Comparison of results on Z(,(10610) and Z(,(10650) parameters obtained 
from T(5S) — > ht(nP)Tv + TY~ (n — 1, 2) analyses. Quoted values are in MeV/c 2 for 
masses, in MeV for widths and in degrees for the relative phase. Relative amplitude is 
defined as a Zb{10(i50) /az b io6io- 



Final state 


/ib(lP)7r + 7r~ 


h b (2P)ir + TT- 


M(Z 6 (10610)) 


10605.1 ±2.2 +-*■;; 


10596 ± 7 + 2 


T(Z b (10610)) 


11 4 +4.5+2.1 
xi - ft -3.9 -1.2 


1fi +16 +13 
1D -10 -4 


Af (Z 6 (10650)) 


10654.5 ±2. 5 + ig 


10651 ±4±2 


r(Z b (10650)) 


20 o+5- 4 + 21 

ZUM -4.7 -5.7 


19 +11 +8 
lz -9 -2 


Rel. amplitude 


1 o+l.O +0.1 
J--S-0.7 -0.5 


-1 ,+3.1 +0.4 
10 -1.1 -0.7 


Rel. phase, 


1SS+ 44 + 4 
loo -58 -9 


9 r,r+56 +12 
iJO -72 -183 



Table 4. Comparison of results on zTj,(10610) and Zf,(10650) parameters obtained 
from T(5S) — > T(nS')7r + 7r~ (n = 1, 2, 3) analyses. Quoted values are in MeV/c 2 for 
masses, in MeV for widths and in degrees for the relative phase. Relative amplitude is 
defined as az b (io650)/az,,io6io- 



Final state 


T(1S)7T + Tr~ 


T(2S)rr + Tr _ 


T(3S)tt + 7t- 


A/(Z b (10610)) 


10609 ± 3 ± 2 


10616 ±2+1 


10608 ± 2 + ^ 


r(Z 6 (10610)) 


22.9 ±7.3 ±2 


21.1 ±4+3 


12.2 ± 1.7 ±4 


M(Z b (10650)) 


10660 ± 6 ± 2 


10653 ± 2 ± 2 


10652 ± 2 ± 2 


r(Z 6 (10650)) 


12 ±10 ±3 


16.4±3.6+J 


10.9 ±2. 6+^ 


Rel. amplitude 


0.59 ± 0.191°;™ 


0.91±0.11±°; O 3 


0.73±O.loi°'o5 


Rel. phase, 


53 ± 6ll 5 50 


-20 ± 18±l 4 


6 ± 2i±f 9 



For the analysis of the T(55) — > T (nS)n + n~ , in addition of the the missing mass 
MM(ir + ir~) associated with the 7r + 7r~ system, the events are further identified by the 
H + H~ pair with an invariant mass in the range of 8.0 GeV/c 2 < Af(/i + /i~) < 11.0GeV/c 2 . 
Figure [16] shows Dalitz plots of the events in the signal regions for the three decay chan- 
nels under study. In all cases, two horizontal bands are evident in the Y(nS)Tr system 
near 10.61 GeV/c 2 (~ 112.6 GeV 2 /c 4 ) and 10.65 GeV/c 2 (~ 113.3 GeV 2 /c 4 ). The 
amplitude analyses of the three-body T(55) — > T(nS)ir + ir~ decays that are reported 
here are performed by means of unbinned maximum likelihood fits to two-dimensional 
Dalitz distributions. The details of the description of the Dalitz plot analysis can be found 
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Figure 16. Dalitz plots for T(nS>+7r- events in the (a) T(1S); (b) T(2S); 
(c) T(3S) signal regions. Dalitz plot regions to the right of the vertical lines are in- 
cluded in the amplitude analysis. 

in 1 93, 94]. Results of the fits to T(55) —$■ T(nS)Tr + ir~ signal events are shown in 
Fig. [17] where one-dimensional projections of the data and fits are compared. To com- 
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bine Z^ and Z^ events we plot T(tiS)tt mass distributions in terms of M(T(nS)w) m i n 
and M(T(n5)7r) max ; fits are performed in terms of M (Y(nS)ir + ) and M(T(nS)n~). 
Results of the fits are summarized in Table |4] The combined statistical significance of the 
two peaks exceeds IOct for all tested models and for all Y(n5)7r + 7r~ decay modes. 





So.4 10.45 10.5 10.55 10.6 10.65 10.7 10.75 



L_£ 




0.58 10.6 10.62 10.64 10.66 10.68 10.7 10.72 10.74 
»{y(3S)Il)_, (CV/c 2 ) 



Figure 17. Comparison of fit results (open histogram) with experimental data (points 
with error bars) for events in the T(15), T(2S) and T(3S I ) signal regions. The hatched 
histogram shows the background component. 
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Figure 18. Comparison of Z;,(10610) and ^(10650) parameters obtained from dif- 
ferent decay channels. The vertical dotted lines indicate B* B and B* B thresholds. 

The Belle experiment observed two charged bottomonium-like resonances, the Zf,(10610) 
and Z;,(10650), with signals in five different decay channels, T(n5)7r ± (n — 1, 2, 3) and 
hb(mP)ir (m = 1, 2). All channels yield consistent results as can be seen in Fig.fTSI 
A simple weighted averages over all five channels give A/[Zf,(10610)] = 10608.4 ± 
2.0MeV/c 2 ,r[Z 6 (10610)] = 15.6±2.5 MeVandM[Z h (10650)] = 10653.2±1.5MeV/c 2 . 
r[Zb(10650)] = 14.4 ± 3.2 MeV, where statistical and systematic errors are added 
in quadrature. The measured masses of these states exceed by only a few MeV/c 2 the 
thresholds for the open beauty channels B*B (10604.6 MeV) and B*B* (10650.2 MeV). 
This "coincidence" can be explained by a molecular-like type of new states, i.e., their 
structure is determined by the strong interaction dynamics of the B*B and B*B meson 
pairs [95]. The T(55) — > hb(mP)ir + tt~ decays seem to be saturated by the Zb(10610) 
and Zb(10650) intermediate states; this decay mechanism may be responsible for the high 
rate of the T(55) — > hb(mP)ir + tt~ process measured recently by the Belle Collabora- 
tion Hi. 



5. Summary 

In this paper, I review the most recent experimental results on the light hadron, charmo- 
nium and bottomonium states. Especially, I gave a review of the exotic candidates ob- 
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served so far. Although a lot of great progresses have been made in the last few years, the 
situation is far from being completely clarified. In most case, new data and high statistics 
are needed, and they are expected to come from the LHC and the future flavor factories. 
At the same time, the analysis effort of the B-Factory experiments is still on going and 
new results are expected to come in the near future. Finally, we need more theoretical 
efforts to make further investigations, for example, in the case of T(55 f ) — > T(nS')7r + 7r _ 
and ht,(nS)Tr + ir~ interpretation and their connections with the e + e~ — > hadrons cross 
sections. In conclusion, hadron spectroscopy is still an intriguing field, new or more pre- 
cise measurements will continue to provide, in the near future, important information to 
better understand QCD and its effective treatments, with a broad impact on many other 
fields. 
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